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Abstract. I survey the use and impact of adaptive mesh refinement (AMR) 
simulations in numerical astrophysics and cosmology. Two basic techniques 
are in use to extend the dynamic range of Eulerian grid simulations in multi- 
dimensions: cell refinement, and patch refinement, otherwise known as block- 
structured adaptive mesh refinement (SAMR). In this survey, no attempt is 
made to assess the relative merits of these two approaches. Rather, the discus- 
sion focuses on how AMR is being used and how AMR is making a scientific 
impact in a diverse set of fields from space physics to the cosmology of the 
early universe. The increased adoption of AMR techniques in the past decade 
is driven in part by the public availability of AMR codes and frameworks. I 
provide a partial list of resources for those interested in learning more about 
AMR simulations. 



1 Introduction 

Since its introduction roughly 20 years ago adaptive mesh refinement 
(AMR) has emerged as an important class of numerical techniques for improv- 
ing the accuracy and dynamic range of grid-based calculations for fluid dy- 
namics problems. Such problems, especially compressible flow, develop steep 
gradients (shock waves and contact discontinuities) which, in the absence of 
mesh refinement, become sources of error for the global solution (e.g., Q). 
Through appropriate local mesh refinement, AMR can be thought of as a 
numerical technique for optimizing the quality of a numerical solution for a 
given computational cost (e.g., 

The last ten years have seen the application of AMR methods to problems 
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and cosmology [ISi llil l2(i l2ll l22L 12,1 . 124 . |2^, '26]. Here, a variety of physical 
processes may operate singly or together in astrophysical fluids to expand the 
range of important length- and time-scales. Such processes include gravity 
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and gravitational instability, reaction kinetics, magnetic reconnection, radia- 
tion transfer, ionization fronts, etc. AMR has also been applied to the solution 
of gravitational N-body problems [27I , and to hybrid particle-fluid simu- 
lations in cosmology 18, 2^. In such applications, the physics is intrinsically 
multi-scale, and AMR can be thought of as a numerical technique for extend- 
ing the dynamic range of resolved physics, regardless of the computational 
cost. It is these sorts of applications which are reviewed here, and where 
AMR can make a scientific and not just an economic impact. 

AMR is also having a positive impact on the methodology of computa- 
tional physics itself. I am referring to the validation of computational codes 
through resolution studies. With AMR, it is now practical on today's super- 
computers to perform resolution studies over a sufficient range of scales to 
obtain convergent results on properties of interest which may be compared 
with laboratory experiments (e.g., 

Table 1 shows the impressive diversity of topics AMR has been applied 
to. Particularly interesting is the range of physical processes AMR has been 
adapted to model. This is listed in the second column of Table 1. As of today, 
AMR has been successfully applied to ideal gas dynamics (Newtonian and 
special relativistic) , reactive gas dynamics, MHD (ideal and resistive), self- 
gravitating gas dynamics and MHD, N-body dynamics, and hybrid fluid/N- 
body systems. As we heard at this conference, methods are under development 
for radiative transfer (Howell), radiation hydrodynamics (Weaver), and solid 
mechanics (Falle) . It is clear AMR is a method of wide applicability, and one 
that is growing in its impact. This growth is fueled in part by the public 
availability of AMR codes and frameworks (see Section 3). 



Topic 


Physics 


Select References 


Code validation 


HD, reactive HD 


[12] 


Solar and space physics 


MHD 




Supernovae and nucleosynthesis 


reactive HD 


[29, 30, 31] 


Interstellar medium 


HD, MHD 


\9, 15, 32, 33] 


Star formation 


grav HD, grav RHD 


fl3, 11] 


Astrophysical jets 


HD, rel HD 


\8, 34, 35] 


N-body dynamics 


particles, grav 


[2^, 2^ 


Hydrodynamic cosmology 


hybrid 


[18, 26j 



Table 1. Classes of AMR applications. 



In this paper I survey the use and impact of adaptive mesh refinement 
simulations in numerical astrophysics and cosmology. Two basic techniques 
are in use to extend the dynamic range of Eulerian grid simulations in 
multi-dimensions: cell refinement (CR), and patch refinement (PR), other- 
wise known as block-structured adaptive mesh refinement (SAMR). Details 
of these two approaches are given elsewhere in this volume. In this review, 
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no attempt is made to assess the relative merits of these two approaches. 
Rather, the discussion focuses on how AMR is being used and how AMR is 
making a scientific impact in a diverse set of fields from space physics to the 
cosmology of the early universe. At the end, I provide a partial list of software 
resources for those interested in learning more about AMR simulations. [NB. 
The following is meant to be representative rather than complete. I apologize 
to authors in advance if their work is not mentioned.] 



2 Applications of AMR 
2.1 Solar and Space Physics 

Coronal mass ejections (CMEs) are transient solar events in which mass and 
magnetic field are ejected from the solar surface. These dynamic events orig- 
inate in closed magnetic field regions of the corona. They produce large-scale 
reconfigurations of the coronal magnetic field and generate solar wind distur- 
bances which manifest as geomagnetic storms here on Earth. Groth et al. 0,0] 
have applied AMR to 3D simulations of solar coronal outflows and ejections. 
The central question which motivates their research is what is the mechanism 
and timescale for coronal mass ejections (CMEs)? The simulations solve the 
equations of ideal MHD in 3D Cartesian geometry, supplemented with the 
Sun's gravitational field and a coronal heating term. The equations are solved 
using the upwind, cell-centered, finite volume scheme of Powell |3fi[l3^ . which 
is a Godunov-type scheme for ideal MHD. Numerical fluxes are computed 
using the approximate Riemann solvers of both Roe and Harten et al. 
|39| . adapted to the MHD eigensystem. The MHD solver is married to the 
parallel AMR framework of de Zeeuw et al. In this approach, the base 
grid is decomposed into blocks of constant size, each of which is assigned to 
a separate processor. If any cell within a block are flagged for refinement, the 
entire block is refined by a factor of two in each dimension, resulting in eight 
sub-blocks of the size of the parent block. 

The CME problem is initialized with a model of the quiescent solar wind 
which involves a high speed wind and open magnetic field lines at high so- 
lar latitudes, and a lower wind speed and closed magnetic field lines in the 
equatorial region. The CME event is triggered with the introduction into the 
computational domain of a density pulse in the closed field region. The mass 
loading inflates the closed field region until it bursts open. Uisng AMR, they 
are able to follow the CME outburst to | AU. According to the authors, the 
benefits of AMR in this demonstration calculation are cost/memory savings, 
and the ability to refine a region of interest (the magnetic X-point) which 
moves through the computational volume. Although they were unable to an- 
swer the ultimate question, the authors argue that AMR combined with more 
realistic solar magnetic field configurations and initiation mechanisms will 
lead to improved understanding of CMEs. 
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2.2 Supernovae and Nucleosynthesis 

An important and growing class of applications for AMR simulations is cal- 
culating the mechanisms and chemical yields of supernovae models in three 
dimensions. In recent years, the importance of convective/turbulent motions 
in both Type la (thermonuclear) and Type II (core collapse) supernovae has 
been recognized 0, \^ . Fully 3D simulations are thus required, and AMR 
serves two fundamental roles. The first is reducing memory and cpu require- 
ments for expensive 3D simulations, and the other is resolving and tracking 
dynamic interfaces where energy generation and nucleosynthesis takes place. 

Type la supernovae are believed to result from the explosive burning of 
carbon and oxygen in a white dwarf (WD) accreting matter from a binary com- 
panion. As the WD's mass approaches the Chandrasekhar limit, a small mass 
increase causes a substantial contraction of the star. The compression raises 
the temperature and thereby thermonuclear reaction rates, liberating energy 
which raises the temperature yet further. Because the WD is degenerate, a 
thermonuclear runaway ensues until thermal pressure becomes comparable to 
the degenerate electron pressure. At this point, the WD begins to expand, 
but not so fast that the thermonuclear reactions are quenched. According to 
models |0, l3]| , a nuclear burning front propagates outward from the center 
of the expanding WD liberating of order 10^^ ergs of energy — sufficient to 
unbind the star. 

Gamezo et al. [s^ have developed a 3D AMR code to simulate the physics 
described above, and have addressed several key questions: (1) How does a 
degenerate C/0 white dwarf explode and with what energy? (2) What is the 
nature and structure of the burning front? And (3) what fraction of the WD is 
burned? They solve Euler's equations of gas dynamics for four reactive species 
coupled to a nuclear reaction network. The energy equation includes nuclear 
energy generation and neutrino losses, as well as electron thermal conduction. 
A key additional incrcdicnt to the physical model is a flame capturing tech- 
nique which adds one additional PDE to be solved for a reaction progress 
variable. The system is evolved for one octant of the WD star on a 3D Carte- 
sian grid using a cell-refinement AMR technique called FTT (Fully Threaded 
Tree) developed by Khokhlov The benefit of cell-refinement over SAMR 
in this application is that the region of interest is a surface — the flame front — 
which is more efficiently captured with a cell-refinement technique il4j . 

Gamezo et al. find that the nuclear flame front is a deflagration front, 
meaning it advances subsonically into the unburned gas. The flame front is 
Rayleigh- Taylor (hereafter, RT) unstable and becomes highly convoluted, re- 
sembling the head of a cauliflower. AMR is used to track this highly distorted 
surface as it advances into the unburned star. They find that a healthy ex- 
plosion results, but that a substantial fraction of the WD remains unburned 
at disassembly. This is contrary to observations , wh ich are found to be more 
consistent with detonating WD models (e.g., |43,|4J|). Gamezo et al. speculate 
that the turbulent flame may trigger a detonation at some point, which would 
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complete the burning and remove this discrepancy. The physics of deflagration 
to detonation transition (DDT) depends on resolving the formation of "hot 
spots" behind the turbulent flame front on scales comparable to the flame 
front thickness This is tiny compared to the radius of the WD, and thus 
not captured by global simulations despite the use of AMR. At present DDT 
must be studied using simulations in small volumes ("subscale simulations".) 

An example of this is the work of Timmes et al. , who used the FLASH 
code ^3 to study the cellular structure of carbon detonation waves in degen- 
erate WDs. Their goal was to assess the effect of numerical resolution on the 
size and shape of the detonation cells and the shock wave interactions that 
create them. FLASH is a block-structured AMR code combining the piecewise 
parabolic method for gas dynamics, nuclear reactions, and the PARAMESH 
AMR library 46]. The simulation was done in a 2D column of gas with con- 
ditions representative of the WD. The detonation was initiated by sending a 
shock wave into the column with strength consistent with a steady ID deto- 
nation wave ( "C J wave" ) . AMR was used to refine by up to a factor of 8 in 
cell size the reaction zone behind the shock. AMR thus provided the ability to 
have high resolution just in the detonation wave and track it as it propagated 
many times its width. 

Timmes et al. were able to converge on the size and shape of the detonation 
cells with rather modest resolution (about 20 cells per burning length scale), 
but found that the size and distribution of pockets of unburned fuel was 
very sensitive to resolution. As long as this scale is small compared to the 
WD radius, observations would average over these compositional differences. 
However, the burning length scale and hence detonation cell size becomes 
comparable to the WD radius as density drops toward the edge of the star ^3 • 
If the spectra of Type la supernovae reflect large scale abundance variations 
due to incomplete combustion, the results of Timmes et al would help define 
the minimum resolution requires to converge on the inhomogeneities. 

In a related work, Zingale et al. have applied FLASH to the dynamics 
of helium detonation on the surface of a neutron star in 2D. In this application, 
as in Gamezo et al., AMR provides the ability to resolve and track the nuclear 
burning front in an expanding envelope of gas. 

AMR is also making an impact in the understanding of the iron core 
collapse supernova explosion mechanism and subsequent explosive nucleosyn- 
thesis. Supernovae classified by observers as Type lb and II are believed to be 
powered by the copious flux of neutrinos emitted as the iron core of a massive 
star implodes to form a neutron star or black hole . Earlier 2D simulations 
showed that neutrino heating sets up convective motions in material accre ting 
onto the proto-neutron star/black hole, breaking spherical symmetry [4lll5fil |. 
The consequence of this is that the shell of 56Ni which forms just outside this 
convective region and is later ejected by the supernovae is highly perturbed 

AMR has not yet been applied to the very difficult problem of neutrino 
transport/heating in a 3D convective flow, although a large collaboration is 
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attempting this [51|. Rather, Kifonidis et al. |29ll52l| have used AMR to sim- 
ulate the evolution of the lumpy nickel shell as it is ejected by the explosion. 
AMR is used to track the expansion of the shell over a factor of 100 in radius 
and its fragmentation by the RT instability while maintaining high resolution 
in the shell. The authors map the results of a 2D axisymmetric core collapse 
simulation 30 ms after bounce into AMRA, a 2D AMR hydro code developed 
by Plewa & Miiller 53] . The equations of multispecies reactive hydrodynamics 
are solved on an adaptive spherical polar grid using the PPM algorithm 
and the Consistent Multifluid Advection (CMA) scheme of Plewa & Miiller. 
The latter minimizes numerical diffusion of nuclear species while preserving 
local mass conservation. Kifonidis et al. find that the 56Ni and other newly 
formed iron group elements are distributed throughout the inner half of the he- 
hum core by RT instabilities operating at the (NH-Si)/0 and (C-|-0)/He shell 
interfaces seeded by perturbations from convective overturn during the early 
stages of the explosion. Using AMR, they are able to carry the calculation out 
to 20,000 sec post-bounce, follow the details of RT growth and mixing. The re- 
sult is that fast-moving clumps of 56Ni are formed with velocities up to 4000 
km/s. This offers a natural explanation for the mixing required in light-curve 
and spectroscopic modeling of Type lb explosions, including SN1987a. 

2.3 Interstellar Medium 

The essential complexity of the interstellar medium (ISM) stems from the fact 
that it is not in equilibrium, dynamically or thcrmodynamically. A variety of 
heating and cooling processes operate to split the ISM into multiple ther- 
mal phases, each with their own characteristic densities, temperatures, and 
evolutionary timescales. Self-gravity concentrates the cold, dense phase into 
molecular clouds, which birth new stars. The most massive stars create large 
amplitude disturbances in the local heating rate through supernova shocks 
and ionization fronts, changing the dynamical and thermodynamical state of 
the gas. To simulate the ISM is rather akin to simulating the Earth's weather, 
which is affected by both local and global influences, and exhibits complexity 
in space and time on a vast range of scales. With AMR simulations, we may 
eventually be able to build an integrated model of the ISM which captures 
both its structural and statistical properties. At present, researchers are look- 
ing at individual processes in greater detail than ever before that ultimately 
may become part of an integrated model. 

The propagation of interstellar shock waves in the ISM has received con- 
siderable attention because of their central role as a source of heat and mo- 
mentum to the interstellar gas. The phenomenology is rich because the ISM 
is inhomogeneous, and the shocks themselves are typically radiative. A set- 
piece calculation is the interaction of a strong planar shock wave striking an 
isolated interstellar cloud, idealized as homogeneous and spherical. The inter- 
action results in the compression and ultimate shredding of the cloud by a 
combination of Richtmyer-Meshkov, Ray leigh- Taylor, and Kelvin-Hclmholtz 
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instabilities. Since these instabilities all grow on the shock-accelerated cloud- 
intercloud interface, it is important that the interface be tracked with high 
resolution. 

The first AMR simulations were performed by Klein, Mc Kee & Colella . 
They wanted to calculate how long a shocked cloud would remain intact before 
mixing into the ISM. They solved the gas dynamical problem in 2D using 
block-structured AMR and a second-order Godunov hydrodynamics scheme 
as described in 0|. They found that for strong incident shocks, the evolution 
is determined by a single parameter: the ratio of cloud to intercloud densities. 
Using AMR, they performed a resolution study to determine the minimum 
resolution needed to capture the most destructive modes of the instabilities. 
They found that a minimum of 100 cells per cloud radius are needed with 
their second order-accurate scheme, and that the cloud is totally fragmented. 

A closely related problem is the long-term evolution of the dense, metal- 
rich clumps ejected by supernova explosions (cf. [l^l). Because dense clumps 
are decelerated less than the diffuse interclump ejecta, they catch up to the 
supernova remnant (SNR) shell and puncture it from behind. The clump first 
encounters the reverse shock, traverses the high pressure intershock region, 
and then exits the forward shock if it survives to encounter the onrushing 
ISM. Cid-Fernandes et al. ^ used the 2D AMR code AMRA [s^ to simulate 
the evolution of a single clump with properties appropriate to a Ni clump 
from a core collapse SN. They initialized the freely expanding envelope of 
a Type II SN in spherical polar coordinates, assuming axisymmetry. They 
placed a clylindrical plug of denser material just upstream of the reverse 
shock, and followed its subsequent evolution by solving the equations of ideal 
gas dynamics including equilibrium radiative cooling. Three levels of mesh 
refinement tracked the blob, providing local resolution equivalent to a uniform 
grid of 1536 x 160 grid cells. The benefit of AMR over a uniform grid was a 
cost savings of 350%. They found that the clump is strongly compressed in the 
intershock region by virtue of the high pressure and strong radiative cooling 
in the clump. While lacking the resolution of Klein et al. |9|, they concluded 
the clump would most likely be disrupted into secondary fragments before 
reaching the dense outer shell. They suggested that x-ray fiares would result 
if the largest of these secondary clumps survived to strike the dense shell. 

Two other AMR-enabled simulations of interstellar shock waves illustrate 
the richness of the phenomena. In the first, Poludnenko, Frank & Blackman 
[ssl simulated the propagation of a planar, adiabatic shock through a clumpy 
medium. The motivation for the simulations was to understand mass-loading 
and mixing of stellar and galactic outflows by inhomogeneities in the am- 
bient medium. A parameter survey was carried out to assess the effects of 
clump mass and spatial distribution on the flow. 2D AMR simulations where 
performed in planar geometry using the AMRCLAW package of Berger & 
LeVeque [s^, which combines SAMR with a second order-accurate Godunov 
scheme for ideal gas dynamics. The contribution of AMR was to achieve high 
resolution in each of many clumps scattered randomly throughout the volume 
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as they are shocked and sheared. Resolution equivalent to a 800 x 1600 uni- 
form grid was achieved in the clumps for a fraction of the cpu and memory 
cost, which is particularly important when conducting a large parameter sur- 
vey. They found that a critical longitudinal and transverse separation between 
clumps exists such that for d < dcrit and L < Lcrit, the post-shock flow is 
strongly interacting, leading to enhanced turbulence and mixing. 

In the second, a ne w ty pe of instability was discovered in radiative shocks 
by Walder and Folini |57| . Radiative shock waves are found in many types 
of classical nebulae, like supernova remnants, planetary nebulae, Wolf-Rayet 
ring nebulae, etc. Typically, a shell of dense material is formed by the inter- 
action of a fast outflow with the circumstellar medium (CSM) or interstellar 
medium (ISM). This shell will be bounded by two shocks: an outer, forward 
shock compressing the CSM/ISM, and an inner, reverse shock compressing 
the ejecta. Typically, the outer shock is radiative, while the inner shock is not 
because of the relative densities involved. The shocked media are separated 
by a contact discontinuity which is normally RT stable at late times because 
the shell decelerates. 

However, it is well known that strongly radiative shocks suffer from an 
overstability such that the shock oscillates about its steady-state position in 
the rest frame of the contact discontinuity ^53|. In multidimensions, Walder 
and Folini showed that different sections of the radiative shock oscillate 
independently, creating lateral pressure disturbances within the high pres- 
sure shell. These chaotic perturbations sometimes accelerate the CD and 
thereby excite the RT instability. This leads to fingers and clumps of dense, 
shocked CSM/ISM efficiently mixing with the shocked ejecta. It is suggested 
by the authors that the mixing will boost the X-ray emission, contribute 
to rapid variability in the emission spectra, and may contribute to the the 
clumpy/filamentary appearance of the nebulae mentioned. The numerical 
simulations were carried out in 2D using AMRCART, which combines block 
SAMR with a second order accurate Godunov solver for gas dynamics. Op- 
tically thin radiative cooling was included assuming equilibrium ionization. 
The simulation was carried out in the rest frame of the CD, and hence the 
role of AMR was not shock tracking, but rather maintaining high resolution 
near the unstable interface. 



2.4 Star Formation 

The gravitational collapse of a gas cloud to form a star is a notoriously difficult 
problem because of the large range of length- and time-scales that need to be 
resolved. This problem has been solved in one dimension assuming spherical 
symmetry using moving adaptive meshes and implicit time integration [sfll 
|S0(. Important structures that need to be resolved span some 9 decades in 
radius from inside the hydrostatic protostellar core to the edge of the accreting 
cloud. Important timescales vary by 6 decades from the sound crossing time 
in the hydrostatic core to the accretion time. AMR holds forth the promise 
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that simulations covering this range of scales will become possible in 3D, 
permitting a self-consistent investigation of disk accretion and the dynamical 
role of magnetic fields. While this is still out of reach, important strides have 
been made on the early st ages of non-spherical cloud collapse by Richard Klein 
and colleagues ■ They are investigating the important problem 

of binary star formation using adaptive mesh techniques. They calculate the 
gravitational fragmentation of slowly rotating molecular cloud cores in 3D, 
assuming the gas is isothermal and non-magnetic. The equations of isothermal, 
self-gravitating gas dynamics is solved using a second order-accurate Godunov 
scheme on a block SAMR grid. The Poisson equation for t he g ravitational 
potential is solved using multigrid relaxation as described in |13|. 

An important early finding was that the cloud would fragment artificially 
due to numerical perturbations unless the local Jeans length is resolved by at 
least 4 cells at all times [l^l- Since the Jeans length scales as p"^^^, if collapse 
raises the central density by a factor a;, then the grid spacing must be decreased 
by a factor in order to avoid artificial fragmentation. The simulations 

of Truelove et al. p^H^ used AMR was to follow compressions of a; 10^ at 
which time the isothermal assumption breaks down. Local mesh refinement of 
up to a factor of 10^ in resolution was accomplished by recursively adding grids 
refined by a factor of 4 wherever the Jeans condition J = Ax/\j < 1/4 was 
about to be violated. Simulations with fixed dynamic range X will inevitably 
violate the Jeans condition when x > X. Truelove ct al. showed that an 
unperturbed cloud will fragment artificially shortly after the Jeans condition 
is violated They went on to show that a slowly rotating cloud with an 
m=2 nonaxisymmetric perturbation does not fragment into a binary system 
as had previously been reported on the basis of fixed resolution simulations by 
Burkert & Bodenheimer |6l|. but rather forms a singular isothermal filament 
in accord with analytic predictions 62] . 

At sufficiently high number densities, however, the isothermal assumption 
breaks down because the cloud becomes optically thick to its own cooling 
radiation [i^- Boss et al. [H^l simulated the non- isothermal evolution of the 
singular filament using both fixed as well as AMR grid codes. Opacity ef- 
fects were modeled by adopting a barotropic equation of state with a variable 
gamma-law. The codes agreed providing the Jeans condition was obeyed in the 
fixed grid run. They found that the filament fragments into a binary system 
whose properties depend sensitively on the equation of state. Since the EOS 
only mocks up radiation trapping in an approximate way, the implication is 
that radiative transfer needs to be modeled self-consistently in future AMR 
simulations of cloud fragmentation. A flux-limited radiation diffusion algo- 
rithm for AMR grids has recently been introduced by Howell & Greenough 
[ssl which is beginning to be applied to star formation [lH . 
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2.5 Astrophysical Jets 

Astrophysical jets are highly coUimated, high speed bipolar outflows powered 
by disk accretion onto compact, gravitating objects. They manifest in a sur- 
prising diversity of systems and length scales, ranging from the pc-long optical 
jets from young stars to the Mpc-long radio jets from active galactic nuclei 
and quasars [63. Regardless of their origin, the jets themselves are interesting 
dynamically and morphologically because they sample both the deep poten- 
tial wells where they were launched, and the interstellar/intergalactic medium 
they propagate through. Because of their high degree of coUimation, the jets 
are believed to be hypersonic, and in the case of extragalactic jets, relativistic. 
A near universal feature of astrophysical jets, whether from stars or galaxies, 
is the occurrence of emission knots. Emission knots are patches of high emis- 
sivity arrayed along the length of the jet. These are generally interpreted as 
internal shock waves in the jet excited by internal or external perturbations 

Beginning in the early 1980s, there have been extensive numerical studies 
of the structure and dynamics of astrophysical jets with the framework of 
ideal gas dynamics and MHD. From the standpoint of simulations, the salient 
difference between protostellar jets and extragalactic jets is that the former 
are dense enough to be strongly radiative, whereas the latter are effectively 
adiabatic. One of the first applications of AMR in astrophysics was by Falle 
& Raga 0, 111 , who studied the detailed structure of an emission knot in a 
radiative protostellar jet. ID models by Raga et al. j62j showed that variations 
in the outflow velocity would create forward-reverse double shock pairs which 
propagate down the jet axis with the mean flow velocity, in good accord 
with observations. However, these models could not capture the bow shock 
appearance of the knots caused by the lateral expansion of the shocked gas 
into the ambient medium. Falle & Raga 8] simulated the multidimensional 
structure of a single knot in the rest frame of the knot. The contribution of 
AMR was to resolve the strong cooling region and ionization structure behind 
the radiative shocks, which can be a small fraction of the jet radius. It is 
impractical to do this with uniform grids. The calculations were performed in 
2D assuming axisymmetry using a block structured AMR grid. Six levels of 
grid refinement were used for an effective grid resolution of 1280 x 640 cells. 
The physics included gas dynamics, solved using the second-order Godunov 
scheme of Falle [7(1 |. non-equilibrium ionization, and radiative cooling. The 
simulations showed that the knots could survive as coherent entities for many 
jet radii, and were morphologically similar to those observed. The ability to 
resolve the ionization structure with AMR allowed them to make synthetic 
emission maps in the commonly detected [S II] doublet and thereby diagnose 
the physical conditions in observed jets. 

Radio observations of radio-loud active galactic nuclei (AGN) mapped 
with VLBI techniques reveal one-sided jets with a stationary core and knots 
of emission that sometimes move superluminally (67| . This is conventially in- 
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terpreted as the result of relativistic Doppler boosting and time dilation when 
observing a relativistic jet at small inclination a ngle s 1711 . With the advent of 
good algorithms for relativistic hydrodynamics \72L l73l l35| , it becomes pos- 
sible to model these sources. The structure, stability, and radio morphology 
of relativistic jets in compact extragalactic radio sources has been studied 
using AMR simulations by Hughes and collaborators 0, |3i Iz^ ■ Given the 
complexities of underlying hydro and the relativistic radiative transfer, they 
simply wanted to know whether such models resemble the data. They car- 
ried out 2D axisymmetric simulations using the code of Duncan and Hughes 
[t^. which combines a second-order Godunov solver for the relativistic gas 
dynamics with the block-structured AMR method of Quirk [t^ . The contri- 
bution of AMR to this work is to crisply capture the internal shocks without 
resort to large uniform grids and supercomputers. Assuming a simple relation 
between synchrotron emissivity and gas pressure, Mioduszuski et al. com- 
puted synthetic radio maps for a variety of inclination angles, Lorentz factors, 
and Mach numbers. They found that the VLBI maps of superluminal sources 
are reasonably well fit with strongly perturbed "pulsed" relativistic jets seen 
nearly end-on. They found that temporal changes of the models' radio ap- 
pearance is not easily related to the underlying hydrodynamic quantities due 
to differential Doppler boosting. 



2.6 Galaxies and Cosmology 

Although AMR was invented for accurately integrating the hyperbolic partial 
differential equations of fluid dynamics, the adaptive mesh used in conjunction 
with the particle-mesh (PM) and P^^M N-body techniques jz^l is extremely 
powerful for simulations of coUisionless, self-gravitating systems of particles as 
arise in galactic dynamics and cosmological structure formation. AMR N-body 
codes have been developed by Couchman Jessop et al. [t^, Kravtsov 
et al. jl^l , Bryan and Norman 0, and Knebe et al. [l^ . These codes 
differ in AMR data structures and how AMR is used to optimize the N-body 
calculation. 

Couchman adaptively introduces one or two levels of block mesh refine- 
ment around highly clustered regions of particles in order to reduce the num- 
ber of particle-particle pairs in a P'^M calculation of dark matter clustering, 
resorting instead to the faster PM calculation on the subgrids. The Poisson 
equation is solved on each level of the grid hierarchy using Fourier techniques 
with special Greens functions. The algorithm is referred to as AP^M (Adap- 
tive P^M). 

Jessop et al. combine the classic PM scheme with block structured local 
mesh refinement to achieve high force resolution in condensed systems. The 
Poisson equation is solved at all levels of the grid hierarchy using Neumann 
boundary conditions interpolated from parent grids and an ADI relaxation 
scheme. The algorithm is referred to as PM^ (Particle Multiple Mesh). 
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Kravtsov et al. combine the octree cell refinement approach of Khokhlov 
Q with PM to create the ART (Adaptive Refinement Tree) code. The Pois- 
son equation is solved using successive multilevel relaxation. 

Bryan & Norman adapted and generalized Couchman's AP'^M to an ar- 
bitrarily deep AMR grid hierarchy and married it to a PPM-derived hydro 
solver (see next section.) Two Poisson solvers are implemented: one based 
on Fourier techniques, and another using multigrid relaxation techniques (see 
O'Shea et al., these proceedings). 

Knebe et al.'s code is similar to the ART code in that cell refinement and 
multigrid relaxation is used, but the underlying data structures and timestep- 
ping schemes are somewhat different. 

Klypin et al. 0| have applied the ART code to the gravitational clustering 
of cold dark matter in cosmological simulations of structure formation. The 
key difficulty of all such calculations is to resolve the scales on which galaxies 
form (1-10 kpc) in cosmological volumes large enough to sample the longest 
perturbation waves or to get good statistics (~ 100 Mpc). The required spa- 
tial dynamic range is therefore 10* — 10^ in 3D for multiple centers of interest 
(galaxies). AMR is one option for achieving such resolution; meshless tree 
codes are another (cf. 80]). Klypin et al. investigated the long-standing over- 
merging problem, in which N-body simulations of the formation of clusters of 
galaxies yield too few galaxy-sized dark matter (DM) halos compared with 
observations Rather, early simulations found that the galaxy DM 

halos merged with one another as they orbited within the cluster potential. 
Originally, it was thought that overmerging was a consequence of the omission 
of dissipative baryons from the models (e.g., 82]). 

Klypin et al. showed that the overmerging problem is primarily a numeri- 
cal resolution problem. Namely, that with inadequate force resolution, galaxy 
DM halos are numerically smeared out. As a consequence, the portion of the 
DM halo that is beyond its tidal radius gets stripped by cluster tidal forces 
as well as through close encouters with other galaxy halos. The DM halos 
essentially evaporate in a few orbits and their cores sink to the center of the 
cluster by dyamical friction. Klypin et al. combined analytic estimates and 
ART simulations to determine the resolution requirements to avoid overmerg- 
ing. The simulations used 128'^ DM particles and a base grid of 256^ cells in 
a volume 15 h""'^ Mpc on a side. Up to 7 levels of 2x cell refinement were 
permitted, for a maximum dynamic range of 32,000 and spatial resolution of 
0.5 h"^ kpc. They found that a force resolution of 1-2 h~^ kpc and a mass 
resolution of ~ 2 x 10^h~^MQ is sufficient to sample the population of galaxy 
DM halos in a rich cluster of galaxies. 

2.7 Hydrodynamic Cosmology 

The marriage of a gravitational N-body code for the cold dark matter in 
the universe with a hydrodynamics code to model the baryonic component 
is referred to as a cosmological hydrodynamics code. A natural marriage is a 
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PM N-body code with an Eulerian gas dynamics code, and a number of such 
codes have been developed jSJ, ISJ, |8a, lla, 113, ■ The spatial resolution of 
these codes is limited to the grid spacing, which limits the spatial dynamic 
range to 1000 or less on current high-end machines. This makes them useful 
for simulations of the diffuse intergalactic medium (e.g., H^), but is far short 
of the 10^~^ dynamic range needed for galaxy large scale structure studies, as 
discussed above. AMR overcomes this limitation. 

AMR hydrodynamic cosmology codes have been developed by Byran & 
Norman 0, 0] , Kravtsov , and Teyssier [oj . The Bryan & Norman code 
Enzo combines a block-SAMR code for ideal gas dynamics using a version of 
the PPM algorithm adapted to cosmological flows [s^, with a PM collision- 
less matter solver as described above (see also paper by O'Shea et al., these 
proceedings.) The code has been supplemented with the multispecies primor- 
dial gas chemistry model of Anninos et al. [o^ , photo- ionization heating by 
an evolving metagalactic UV background, and a parameterized model for star 
formation and feedback. The Kravtsov code builds upon the ART N-body 
code described above, and adds the second-order Godunov solver for ideal gas 
dynamics described in The RAMSES code j^^l developed by Teyssier is 
similar to the Kravtsov code with minor differences in implementation. 

A problem all three groups have attacked is the formation of an X-ray 
cluster of galaxies, treating the baryons as non-radiative. This is a good ap- 
proximation since the 10* K gas characteristic of many X-ray clusters has a 
cooling time long compared to the Hubble time. A particular set of initial con- 
ditions known as the Santa Barbara cluster has served as a community test 
problem, and is described in Frenk et al. The chief difficulty is resolving 
the X-ray core radius (~ 100 kpc) in the forming cluster within a simulation 
volume 64 Mpc on a side. Most of the X-ray luminosity is contained within 
this region. If one wants to resolve the core radius with 10 cells, say, then a 
dynamic range of 6,400 is required. Bryan & Norman |l8l | achieved a dynamic 
range of 8,192 using a base grid of 128'^ cells and 6 levels of 2x refinement. 
These results were compiled with those of 11 other codes and presented in 
[q^ . It was found that to compute the X-ray luminosity to within a factor of 
2 accuracy, at least this resolution is required. Kravtsov, Klypin & Hoffman 
[q^ I simulated the Santa Barbara cluster with the ART code with the same 
resolution as Bryan & Norman, and found excellent agreement with their re- 
sults. These simulations, as well as those of Teyssier |91|, have shown that the 
distribution of thermal pressure in the cluster is a much more robust quan- 
tity. The significance of this is that Sunyaev-Zeldovich (SZ) effect in clusters of 
galaxies, which is proportional to the line-of-sight integral of the intracluster 
gas pressure, is robustly predicted. This makes AMR simulations a powerful 
tool for guiding upcoming observational surveys of high redshift clusters using 
the SZ effect (e.g., 113). 

A second fruitful application of AMR cosmological hydrodynamics con- 
cerns the formation of the first bound objects and stars in the universe. Within 
the CDM model of structure formation, dark matter begins clustering on small 
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mass scales after the epoch of matter-radiation equahty — about 30,000 years 
after the Big Bang. The characteristic mass scale for DM halos increases with 
time such that by redshifts of z=20-30, it becomes comparable to the Jeans 
mass in the expanding, adiabatically cooling, primordial gas. Abel, Bryan & 
Norman have used AMR to simulate how baryons collect into the poten- 
tial well of such a low mass DM halo and the ensuing cooling and contraction 
of the gas to form a primordial molecular cloud in the halo's center. The simu- 
lation used a 64'^ base grid and 12 levels of 2x refinement for a dynamic range 
of 2.6 X 10^. In addition to dark matter, gravity and gas dynamics, the cal- 
culation solved a 9-species chemical reaction network to model the gas phase 
reactions which produce molecular hydrogen — the primary coolant in primor- 
dial gas. At this resolution, a primordial molecular cloud of size 5 pc was 
well resolved in a simulation volume 128 kpc (comoving) on a side at z=19. 
At the end of the calculation, a single, gravitationally unstable cloud core of 
mass ^ IOOM0 began collapsing. To follow the evolution of the collapsing 
core to higher densities, Abel, Bryan & Norman used an additional 15 
levels, for a dynamic range of 10^". The mesh refinement was driven by the 
Jeans condition and an analogous condition based on the local cooling time. 
At this resolution, they were able to refute the prediction by Silk 96] that a 
chemo-thermal instability would fragment the collapsing core into low mass 
stars. At the end of the calculation, only a single, collapsing, fully molecular 
cloud core was found with a size comparable to the Solar System. With mean 
density of 10^^ cm~'^, this core would trap its cooling radiation and become 
hydrostatic. Based on the accretion rate at the end of the simulation, it is 
predicted that the cloud envelope of would accrete in 10"* years, forming a 
Population III star with a mass in the range 30 — 300Mq. A calculation with 
34 levels of refinement (dynamic range 10^^) by Bryan, Abel & Norman |9^ 
confirms this result. 



Code 


Description 


URL 


AMRCLAW 


SAMR infrastructure and hyperbolic solvers 


[98] 


AMRCART 


SAMR application: 3D MHD 


[99] 


BEARCLAW 


SAMR infrastructure and PDE solvers 


flOO] 


CHOMBO 


SAMR infrastructure and PDE solvers 


flOl] 


Enzo 


SAMR application: hydrodynamic cosmology 


fl02] 


FLASH 


SAMR application: reactive fluid dynamics 


fl03] 


MLAPM 


AMR application: cosmological N-body 


fl04] 


NIRVANA 


SAMR application: 2D and 3D MHD 


[105] 


PARAMESH 


SAMR infrastructure 


fl06] 


SAMRAI 


SAMR infrastructure 


fl07] 



Table 2. Downloadable AMR Software. 



3 AMR Software 



Impact of AMR 15 



Here I tabulate some AMR libraries and application codes that are available 
for download (Table 2). This list is incomplete, because of the lack of central- 
ized information about such tools, as well as the rapid rate of development in 
the field. 
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